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ABSTRACT 


Algorithms  for  blind  equalization  and  data  recovery  of  orthogonal  frequency-division 
multiplexed  (OFDM)  signals  transmitted  through  fading  channels  are  implemented  and 
simulated  in  this  thesis.  The  channel  is  estimated  without  knowledge  of  the  transmitted 
sequence  (i.e.,  blindly)  using  a  least  mean  squares  (LMS)  adaptive  fdter  and  filter  bank 
precoders.  This  method  was  used  to  estimate  channel  characteristics  using  both  binary 
and  quadrature  phase-shift  keying  signals.  Additionally,  the  method  was  analyzed  for 
robustness  with  a  poor  initial  estimate  of  channel  characteristics,  with  the  addition  of 
white  Gaussian  noise  to  the  signal,  and  with  non-stationary  channel  conditions. 

In  addition,  it  was  shown  that  the  proposed  method  is  particularly  suited  in 
situations  with  deep  fading  channels,  where  some  of  the  subcarriers  have  a  very  low 
signal-to-noise  ratio  (SNR). 

Simulations  for  both  aspects  of  this  thesis  were  conducted  using  MATLAB,  and 
the  results  are  presented. 
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EXECUTIVE  SUMMARY 


Orthogonal  frequency-division  multiplexing  (OFDM)  modulation  has  become  the 
preferred  modulation  method  for  numerous  modem  standards,  from  IEEE  802.11a/g/n 
(Wi-Fi)  to  IEEE  802.16  (WiMAX)  to  4G  Long  Tenn  Evolution  [1],  The  reason  that 
OFDM  modulation  has  become  so  prevalent  is  because  of  a  number  of  advantages  it  has 
over  other  modulation  techniques. 

One  of  the  main  advantages  of  OFDM  modulation  is  the  ease  of  channel 
equalization,  especially  in  multipath  environments  [1].  This  is  achieved  by  dividing  the 
channel  bandwidth  into  narrow  sub-channels  with  orthogonal  subcarriers,  and 
transmitting  a  M  -ary  quadrature  amplitude  modulation  (MQAM)  symbol  on  each  of  the 
subcarriers.  Even  in  the  presence  of  severe  multipath,  each  subchannel  is  affected  by 
flat-fading,  which  can  be  easily  equalized.  A  guard  time  interval  between  symbols 
prevents  inter-symbol  interference.  Additionally,  OFDM  modulation  supports  high  peak 
data  rates  and  non-line-of-sight  applications  [1]. 

As  OFDM  symbols  are  transmitted  through  a  channel,  they  undergo  distortion  due 
to  noise  and  multipath  effects.  To  counteract  the  effects  of  channel  distortion,  an 
estimate  of  channel  characteristics  must  be  calculated.  In  most  standards,  this  is  done 
using  preambles  and  pilot  signals,  resulting  in  a  loss  of  data  rate.  Additionally,  most 
standards  use  a  cyclic  prefix  (CP)  during  their  guard  intervals,  which  results  in  a  loss  of 
power  efficiency  since  this  information  is  disregarded  at  the  receiver. 

A  method  of  blind  channel  equalization  that  can  be  applied  to  OFDM  has  been 
proposed  by  P.  P.  Viadyanathan  et  al  [2],  [3].  This  method  uses  a  zero  prefix  (ZP)  vice  a 
CP,  the  reason  being  that  the  transmitted  ZP  constrains  the  received  signal  to  a  subspace 
from  which  the  channel  information  is  extracted.  This  increases  both  data  throughput  and 
power  efficiency. 

The  method  of  blind  equalization  proposed  in  [2]  is  investigated  in  this  thesis. 
The  method  is  implemented  in  MATLAB  with  both  off-line  and  real-time  channel 
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estimation  investigated.  Although  the  method  in  [2]  uses  a  ZP,  we  show  that  it  can  be 
extended  to  the  most  traditional  CP  based  OFDM  by  some  simple  preprocessing  of  the 
received  data. 

The  method  proposed  in  [2]  was  implemented  as  a  recursive  algorithm  using  a 
least  mean  squares  (LMS)  adaptive  estimation  method.  However,  it  has  been  found  that 
it  is  very  sensitive  to  the  SNR  of  the  received  signal. 

In  this  research,  the  approach  in  [2]  has  been  re-derived  and  a  third,  alternative 
adaptive  algorithm  implemented.  This  method  also  uses  a  LMS  adaptive  estimation 
method;  however,  it  is  shown  to  be  more  robust.  Both  adaptive  methods  can  be  fine- 
tuned  by  properly  setting  their  step  size.  This  thesis  investigated  two  LMS  step  sizes, 
determined  empirically. 

The  performance  of  these  algorithms  under  different  conditions  is  also  compared. 
The  first  simulation  was  a  comparison  of  their  performance  in  a  channel  with  additive 
white  Gaussian  noise  (AWGN).  The  third  method  of  channel  estimation  with  the  larger 
of  the  two  step  sizes  performed  the  best  in  high  signal-to-noise  ratio  (SNR)  conditions 
(around  30dB).  The  second  simulation  used  a  non-stationary  channel  without  AWGN 
modeled  as  a  simple  random  walk.  In  this  case,  the  original  method  of  channel 
estimation  and  the  alternative  adaptive  method  had  similar  results,  although  the 
alternative  adaptive  method  with  the  larger  step  size  performed  slightly  better  at  channel 
rates  of  change  of  0.08  percent  of  the  total  channel  power. 

Since  both  adaptive  methods  require  an  estimate  of  the  channel  as  an  initial 
condition,  the  final  test  run  on  the  two  algorithms  was  a  comparison  of  how  long  it  would 
take  the  channel  estimate  to  converge  with  an  inaccurate  initial  channel  estimate.  In  this 
case,  the  more  robust  alternative  adaptive  method  with  the  larger  step  size  once  again 
outperformed  the  initial  adaptive  method. 

Finally,  this  thesis  investigates  a  method  of  estimating  data  lost  in  nulls  caused  by 
deep  fading  in  the  wireless  channel,  due  to  multipath.  The  algorithm  was  implemented  in 
MATLAB,  and  compared  to  a  normal  OFDM  receiver  that  does  not  conduct  null 
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estimation.  The  proposed  null-estimating  algorithm  has  similar  error  rates  to  the  original 
receiver  for  SNRs  less  than  approximately  21  dB  but  then  significantly  outperforms  it 
and  recovers  the  data  lost  in  nulls  for  larger  SNRs. 

[1]  J.  G.  Andrews  et  ah,  Fundamentals  ofWiMAX.  San  Francisco,  CA:  Prentice  Hall, 
2007. 

[2]  P.  P.  Vaidyanathan  and  B.  Vrcelj,  “A  frequency  domain  approach  for  blind 
identification  with  filter  bank  precoders,”  in  Proc.  2004  Int.  Symp.  Circuits  and 
Systems,  pp.  Ill-  349-52,  2004. 

[3]  P.  P.  Vaidyanathan  and  B.  Su,  “Remarks  on  certain  new  methods  for  blind 
identification  of  FIR  channels,”  in  Conf  Rec.  38th  Asilomar  Conf.  Signals, 
Systems  and  Computers,  Asilomar,  CA,  pp.  I-  832-836,  2004. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Orthogonal  frequency-division  multiplexing  (OFDM)  modulation  has  become  the 
modulation  method  of  choice  for  wireless  broadband  communications.  It  is  currently  in 
use  in  numerous  standards,  including  IEEE  standards  802.11  and  802.16,  and  military 
wireless  communications  such  as  the  Wideband  Networking  Waveform  of  the  Joint 
Tactical  Radio  System.  The  reason  for  its  proliferation  comes  from  the  advantages  that  it 
has  over  other  methods  of  modulation.  The  main  advantage  is  that  the  technique  allows 
high  data  rate  transmissions  to  be  sent  over  channels  that  experience  multipath  fading. 

One  of  the  challenges  of  this  modulation  technique  is  that  it  requires  knowledge 
of  the  multipath  channel  to  perfonn  equalization.  How  to  acquire  knowledge  of  the 
channel  using  blind  equalization  is  examined  in  this  thesis.  This  is  significant  since  it 
increases  the  throughput  of  the  channel  by  reducing  (or  eliminating)  the  need  for  pilot 
signals  normally  embedded  in  the  OFDM  symbol. 

B.  OBJECTIVES  AND  METHODOLOGY 

The  goal  of  this  thesis  was  to  further  investigate  the  method  of  blind  channel 
equalization  proposed  by  P.  P.  Vaidyanathan,  et  al.  [1].  Once  the  method  had  been 
implemented  in  code,  a  binary  or  quadrature  phase-shift  keyed  (BPSK  or  QPSK 
respectively)  signal  was  sent  through  either  a  channel  with  additive  white  Gaussian  noise 
(AWGN)  or  a  non-stationary  channel  modeled  by  a  simple  random  walk. 

Once  this  was  completed,  the  method  was  implemented  using  a  least  mean 
squares  (LMS)  estimation  technique.  This  method  was  also  tested  using  BPSK  and 
QPSK  signals  through  either  a  channel  with  AWGN  or  a  non-stationary  channel.  It  was 
also  tested  using  an  inaccurate  initial  channel  estimate  to  determine  how  long  the  channel 
estimate  would  take  to  converge  to  the  correct  channel. 
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Additionally,  a  secondary  LMS  method  was  implemented  to  overcome  ill- 
conditioning  issues  with  the  initial  LMS  method.  This  was  tested  in  the  same  manner  as 
the  first  LMS  method. 

Finally,  this  paper  shows  an  additional  advantage  of  the  secondary  LMS  method; 
namely,  that  this  technique  allows  for  the  recovery  of  information  sent  through  deep 
fading  channels. 

C.  BENEFITS  OF  STUDY 

Demand  for  data,  specifically  multimedia  content,  continues  to  increase  in  a 
nearly  exponential  manner  [2].  This  has  led  to  an  increased  demand  on  mobile 
broadband  technologies,  including  4G  phone  technology  and  WiMAX  [3].  Both  these 
emerging  technologies  use  OFDM  modulation.  By  reducing  or  eliminating  pilot  signals 
in  the  OFDM  symbol,  the  amount  of  data  throughput  is  increased.  Additionally,  if  a  zero 
prefix  (ZP)  is  used  instead  of  a  cyclic  prefix  (CP),  it  saves  power  at  the  transmitting 
station,  whether  it  is  a  mobile  or  base  station  [3]. 

D.  ORGANIZATION 

This  thesis  is  organized  into  six  chapters.  Chapter  II  is  a  background  discussion 
of  OFDM  modulation,  and  examples  of  its  use  in  current  standards  are  provided.  Chapter 
III  contains  an  in-depth  discussion  of  the  blind  equalization  method  proposed  in  [1],  An 
adaptive  estimation  method  based  on  the  discrete  Fourier  transform  (DFT)  is  proposed  in 
Chapter  IV.  Simulation  results  from  the  adaptive  method  are  compared  with  simulation 
results  from  the  original  in  Chapter  IV.  A  method  of  recovering  the  information  lost  in 
deep  fading  channels  is  discussed  in  Chapter  V,  and  the  results  of  the  implementation  of 
this  method  are  presented.  Finally,  a  summary  of  the  work  in  this  thesis,  significant 
results,  and  recommendations  for  future  work  are  given  in  Chapter  VI. 
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II.  OVERVIEW  OF  OFDM 


Orthogonal  frequency-division  multiplexing  is  a  method  of  multicarrier 
modulation  used  in  a  variety  of  standards  for  both  wired  and  wireless  communications. 
A  brief  overview  of  OFDM  modulation  and  its  use  in  IEEE  standards  802.11a,  802.1  lg, 
and  802. 1  In  is  provided  in  this  chapter. 

The  principle  behind  OFDM  is  the  idea  that  a  wideband  frequency  channel  can  be 
divided  into  subchannels  with  narrower  bands  on  which  data  can  be  transmitted. 
Orthogonality  of  the  subcarriers,  in  a  sense  specified  later,  prevents  interference  between 
the  subbands.  This  technique  is  particularly  useful  when  transmitting  through  frequency- 
selective  channels  arising  from  environments  that  have  severe  multipath,  and  it  allows 
simple  channel  equalization  at  the  receiver  [4], 

In  OFDM,  the  data  is  transmitted  in  blocks  via  a  sequence  of  OFDM  symbols, 
separated  by  a  guard  time  interval  to  prevent  inter-symbol  interference  (ISI).  This  is 
guaranteed  by  choosing  the  guard  time  between  OFDM  symbols  to  be  larger  than  the 
maximum  channel  time  spread  [4]. 

A  key  issue  addressed  in  this  thesis  is  what  is  best  to  transmit  during  the  guard 
time.  In  the  most  common  implementations  of  OFDM,  such  as  in  IEEE  Standards 
802.1 1,  802.16,  and  others,  the  CP  is  used,  which  is  a  cyclic  repetition  transmitted  as  the 
guard  interval.  A  consequence  of  choosing  to  transmit  the  CP  is  that  it  becomes  very 
easy  to  demodulate  the  signal  in  the  frequency  domain,  since  linear  convolution  becomes 
circular  convolution  where  the  fast  Fourier  transfonn  (FFT)  plays  a  dominant  role. 

In  a  number  of  papers  [1][5],  it  is  shown  that  a  ZP,  obtained  by  transmitting  all 
zeros  during  the  guard  interval,  is  an  attractive  alternative  to  a  CP  because  it  yields 
sufficient  information  to  perform  blind  equalization.  In  other  words,  in  principle  at  least, 
with  the  ZP  the  channel  can  be  estimated  from  the  received  data  without  the  need  of 
pilots  and/or  training  signals.  This  is  due  to  the  fact  that,  with  the  addition  of  the 
transmitted  zeros,  the  received  signal  is  constrained  on  a  subspace  that  yields  the  channel 
information,  regardless  of  the  transmitted  data. 
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This  research  takes  this  issue  even  further.  It  shows  that  another  advantage  of 
using  an  OFDM  symbol  with  ZP  (rather  than  the  standard  CP)  is  the  fact  that  by  using  a 
longer  FFT  at  the  receiver,  subcarriers  with  low  signal-to-noise  ratio  (SNR)  due,  for 
example,  to  deep  fading,  can  still  be  reliably  recovered  using  simple  properties  of  the 
FFT.  Of  course,  this  comes  with  some  added  complexity  but  is  still  manageable  with 
today’s  technology. 

Interestingly,  with  some  minor  preprocessing,  the  ZP  OFDM  approach  can  be 
extended  to  the  existing  CP  OFDM  with  the  same  advantages  of  blind  equalization  and 
robust  recovery  in  deep  fading  channels. 


A.  OFDM  SYMBOL 

In  OFDM,  the  transmitted  signal  is  defined  as 

s(t )  =  Re  |  ej2yTFc‘x(t)^ 

where  Fc  is  the  carrier  frequency,  and  the  baseband  complex  signal  is  defined  as 


(1) 


x(t)=  Y,  ckejl7rkm  0<t<Ts 


symbol 


(2) 


k=— 
k*  0 


where  the  coefficients  ck  depend  on  the  transmitted  data  and  N F  is  the  number  of 
subcarriers.  The  symbol  length  Tnmbol  is  made  of  the  data  length  Th  and  the  time  guard 
length  T  as 


T  —T  +  T 

1  symbol  1  b  '  1  g  * 


(3) 


To  ensure  that  the  subcarriers  are  orthogonal  to  one  another,  the  difference 
between  subcarrier  frequencies  is  chosen  to  be 


A  F  =  —  . 


(4) 


This  results  in  the  following  integral,  which  defines  the  orthogonality  of  the 
subcarriers: 


—  J  ej2KFkt e^j27lF,t dt  =  —  J  eilK<k-l)mdt  = 


1 


Lb 


Lb  t„ 


[1  ifk=l 
lo  if  k  ib  l 


(5) 
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If  the  channel  is  linear  time-invariant  (LTI),  at  least  within  the  symbol  duration, 
and  the  guard  time  Tg  is  greater  than  the  channel  time  spread,  the  subcarriers  are  still 

orthogonal  at  the  receiver  and  are  given  by 

Tg<t<  Tsymit)l ,  (6) 

k 

where  H(F)  is  the  frequency  response  of  the  channel  [4], 


The  discrete  time  implementation  is  very  straightforward.  Let  Fs  be  the  sampling 
frequency,  N  be  the  number  of  data  samples  in  each  symbol,  and  A F  =  1/  AT  =FjN  be 
the  subcarrier  spacing.  Then  (2)  becomes 


Ad 

£ 

N 


jlnk^n-L)  l 


n]% 


jAn-L) 
o  N 


n  =  0.....L  +  N -l 


(7) 


where  the  guard  time  is  Tg  = LTS  and  the  scale  factor  of  1  / N  is  immaterial  but  is  included 

for  later  convenience.  This  result  leads  itself  to  using  the  inverse  fast  Fourier  transform 
(IFFT)  to  complete  the  precoding  necessary  for  OFDM: 


1  2 

x[n  +  L]  =  —  cke 

N  k=-& 

2 


2x 

ik - n 

N 


(8) 


If  X[k]  =  ck  for  positive  subcarriers  (k>  0),  X[N  +  k]  =  ck  for  negative 
subcarriers  ( k  <  0  ),  and  A[/t]  =  0  for  k  -  0 ,  then  (8)  can  be  written  in  terms  of  the  IFFT 
as 


i  N- 1  y  2;r 

x[n  +  L]  =  -YJX[k]e  N"  =IFFT{X[k]}.  (9) 

X  k= o 

From  (2)  and  the  subcarrier  frequencies  being  multiples  of  A F  in  (4),  it  is  easy  to 
see  that  the  data  transmitted  during  the  guard  interval  is  just  a  periodic  repetition.  In 
particular,  the  CP  is  comprised  of  the  last  L  points  of  the  IFFT  { A[A]| . 


The  demodulation  of  the  OFDM  signal  is  based  on  the  convolution  of  the  channel 
impulse  response  and  the  signal: 
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(10) 


1. 


y[n ]  =  h[n]*x[n] 

i  JV-l 

=  h[n]*~YJX[k]e 

-/V  k=[ 


.2  7T 

i - kn 

N 


k= 0 


=  77  Y,H\k]X[k]e 

N  k= 0 

=  IFFT  {H[k]X[k]} 


,2n 

i - kn 

N 


Based  on  the  results  of  (10),  the  following  conclusion  can  be  drawn: 
H[k]X[k\  =  FFT  { y[L\,  ...,y[L  +  N- 1]} 
with  H[k]  =  FFT  {//[0], . . . ,  h[L  - 1],  0, . . . ,  0}  it  =  0,...,  JV-l’ 


(11) 


The  block  diagram  of  a  simple  OFDM  communication  system  is  shown  in  Figure 


r  *10] 

All) 

AfiV-ii 


r  WlO]A[0] 
W[1)A'[11 

: 


Figure  1.  A  block  diagram  of  a  simple  OFDM  communications  system.  From  [4], 

Since  for  the  k  -th  subcarrier,  the  received  signal  is  of  the  form 

Y[k\  =  H[k\X[k\  +  W[k'\  (12) 

where  W[k]  is  the  FFT  of  vtf/z] ,  additive  noise,  the  transmitted  signal  X[k]  can  be 
recovered  using  a  Weiner  Filter  as 
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(13) 


X[k]  = 


H*\k\ 


\H[k]\ 


Y\k] 


+  cr: 


if  the  channel  frequency  response  is  known. 


There  are  some  limitations  to  this  approach  of  OFDM  modulation.  The  CP  does 
not  contain  useful  data  and  is  just  overhead  for  the  modulation  technique.  Additionally, 
the  channel  imposes  some  constraints  on  the  length  of  the  data  interval  such  that  it  must 
be  larger  than  the  maximum  channel  time  spread  to  minimize  overhead  and  must  be 
much  smaller  than  then  period  associated  with  the  Doppler  spread  to  maintain  the  LTI 
channel  assumption  [4]. 


The  advantages  of  OFDM  are  that  it  is  effective  in  multipath  environments,  it  is 
efficient  in  implementation,  and  the  data  rate  can  be  optimized  based  on  channel 
conditions  and  noise  characteristics.  The  disadvantages  are  that  it  is  sensitive  to  time  and 
frequency  synchronization,  and  data  can  be  lost  when  there  are  spectral  nulls  in  the 
channel  frequency  response  [4]. 


B.  OFDM  IN  STANDARDS 

1.  IEEE  Standard  802.11a 

IEEE  Standard  802.11a,  which  describes  the  high-speed  physical  layer  in  the 
5GHz  band,  uses  OFDM  to  transmit  data  over  wireless  local  area  networks  (LANs)  [6]. 
The  standard  uses  52  sub-carriers,  four  of  which  carry  the  pilot  signal  and  48  that  carry 
the  data.  The  A F  is  set  at  0.3125  MHz,  with  the  guard  interval  0.8  /is  long  and  the  data 
interval  3.2  jus  long.  The  mapping  of  the  52  sub-carriers  to  the  time  domain  outputs  can 
be  seen  in  Figure  2. 
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Tune  Domain  Outputs 


Figure  2.  Mapping  of  IFFT  inputs  to  time-domain  outputs  for  802.1  la.  From  [7]. 

2.  IEEE  Standard  802.1  lg 

IEEE  Standard  802.1  lg,  is  a  modification  of  the  802.11  standards  for  the  2.4-GHz 
band  to  include  OFDM  modulation.  The  inputs  to  the  IFFT  retain  the  same  structure  as 
those  in  802. 1  la  but  with  a  lower  carrier  frequency 

3.  IEEE  Standard  802.11n 

IEEE  Standard  802.1  In  uses  OFDM  modulation  in  a  multiple-input  multiple- 
output  (MIMO)  scheme  to  attain  higher  throughput  and  more  robustness  to  multipath  and 
noise  for  both  the  2.4-  and  5-GHz  bands  [8].  It  provides  up  to  four  spatial  streams  of 
OFDM  modulated  data  to  achieve  a  maximum  data  rate  of  600  Mbps. 

4.  Other  Standards 

OFDM  modulation  is  also  used  in  various  other  applications.  WiMAX  is  based 
on  IEEE  Standard  802.16,  an  OFDM  based  system.  Other  standards  that  use  OFDM 
modulation  include  power-line  communications,  various  digital  television  and  radio 
standards,  and  the  Wideband  Networking  Waveform  of  the  Joint  Tactical  Radio  System. 
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These  applications  of  OFDM  use  pilot  signals  and  training  symbols  to  conduct 
channel  equalization.  In  the  next  chapter,  we  see  how  channel  equalization  can  be 
conducted  without  the  use  of  these  methods. 
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III.  BLIND  CHANNEL  ESTIMATION  IN  OFDM  WITH  ZERO 

PREFIX 


A  class  of  blind  equalization  algorithms  has  been  presented  by  P.  P.  Vidyanathan 
et  al.  in  a  number  of  papers  [1],  [5],  which  can  easily  be  extended  to  OFDM  with  ZP. 
Their  argument  is  presented  in  this  chapter.  An  alternative  derivation  of  their  results  is  at 
the  basis  of  further  properties  of  this  algorithm,  as  is  presented  in  Chapters  IV  and  V. 

A.  BLIND  EQUALIZATION  USING  A  ZERO  PREFIX 

In  the  method  described  in  [1],  the  CP  used  in  the  standards  discussed  in  Chapter 
II  is  replaced  with  the  use  of  a  ZP.  Instead  of  filling  the  guard  interval  with  data  from  the 
end  of  the  data  block  or  the  last  L  samples  of  the  data  block,  the  data  is  zero-padded  with 
L  zeros  during  the  guard  interval.  As  long  as  the  guard  interval  is  still  longer  than  the 
maximum  time  spread  of  the  channel,  there  is  no  1ST 

In  ZP  OFDM  the  data  is  transmitted  as  a  sequence  of  discrete  time  OFDM 
symbols  x[0],x[l defined  as 

x0[n] 

Xi  [n] 

x[n\=  xMJn\  e  C'’Xi  (14) 

0 

0 

where  M  is  the  length  of  the  FFT  and  IFFT,  and  P  -  M  +  L  is  the  overall  length  of  the 
OFDM  symbol,  including  the  ZP.  Then  each  input  block  is  in  a  subspace  of  Cftl .  If 
there  is  no  noise,  the  received  signal 

ToM 

TiM 

JV-iM 


li 


is  also  in  a  subspace  of  Cftl .  With  this  in  mind,  it  is  possible  to  show  that  the  channel 
impulse  response  c[/]  can  be  estimated  from  the  subspace  of  y[n]  without  knowledge  of 

the  input  signal  x[»] . 

B.  CHANNEL  ESTIMATION 

The  general  principles  behind  the  method  proposed  in  [1]  are  discussed  in  this 
section.  Let  the  n  -th  data  block  be  defined  as 

s0[n] 

s[n]=  i  .  (16) 

_SM- M_ 

The  data  block  is  then  encoded  with  a  matrix  Rt  e  CMxM  which  in  the  cases  of 
OFDM  modulation  is  associated  with  the  IFFT. 

s[n] 

Figure  3.  The  signal  s[«]  is  encoded  with  the  matrix,  then  L  zeros  are  appended. 

After  [7]. 

After  precoding,  the  transmitted  blocks  x[«]  are  defined  as: 


(17) 

(18) 
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with  the  zeros  indicating  the  zero  prefix.  Although  the  appended  zeros  really  represent  a 
suffix,  for  clarity  of  discussion,  they  will  continue  to  be  referred  to  as  the  zero  prefix. 

In  the  section  that  follows,  it  can  be  seen  that  the  received  signal  has  a  particular 
structure  from  which  the  channel’s  characteristics  can  be  extracted  without  knowledge  of 
the  transmitted  data. 


1.  Received  OFDM  Symbol 


The  output  from  the  channel  in  the  time  domain  is  the  convolution  of  the  impulse 
response  of  the  channel  with  x(t) .  Since  there  is  no  ISI  due  to  the  ZP,  each  block  can  be 
examined  separately,  and  the  received  signal  for  the  n  -th  block  becomes 

yi[n\  =  YJc[k]xl_k[n]  l  =  0,1,...,  P-1  (19) 

k=  0 

with  c[k],  k  =  0, . . . ,  L  - 1 ,  the  impulse  response  of  the  channel.  Then  it  is  easy  to  verify 


that 


'  c[0] 

0 

0 

c[l] 

c[0]  ••• 

ToM 

c[l]  -• 

TiM 

c[L- 1] 

0 

0 

c[L~  1]  - 

c[0] 

yP-M 

0 

c[l] 

0 

0 

c[L- 1] 

x0[n] 
xm- iM 


which  can  be  represented  in  matrix  from  as 

y[n]  =  Axd[n]  =  ARxs]ji\ 


(20) 


(21) 


with  y[n]  a  Px  1  vector,  xd[n~\  ,  the  data  only  version  of  *[»],  a  Mx  1  vector,  and  A  a 


PxM  matrix.  The  matrix  A  contains  the  information  of  the  channel  and  it  has  an 
interesting  property.  Given  any  complex  variable  z  ,  it  can  be  shown  that 
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\  z-1. 


z-2, 


'  c[ 0] 

0 

0 

fill 

c[0]  ••• 

fid  ••• 

c[L~  1] 

0 

0 

c\L- 1]  ••• 

c[0] 

0 

c[l] 

0 

0 

c[L~  1] 

=  [C(z),  z-'C(z),  z~2C(z),  z-("-»C(z)] 


,(22) 


with  C(z) ,  the  transfer  function  of  the  channel,  defined  as 


L- 1 


C(z)  =  Xc[/]z-'=Z{c[/]}. 


(23) 


1=0 


The  right  hand  side  of  (22)  becomes 

C(z)[  1,  z-1,  ...,  z-^].  (24) 

Using  the  results  above,  we  see  how  the  transmitted  data  y[»]  can  be  recovered 
from  the  received  data  y[n\  assuming  that  there  is  no  noise.  Let 

s[n\  =  Ey[n\,  Ee  CMxP ,  (25) 

with  E ,  an  arbitrary  M  xP  matrix  to  be  determined  later.  Substitute  (21)  into  (25),  and 
the  following  equation  results: 

y[n]  =  EARvs[n\ 

r  i 

(26) 


=  [  E  ] 


A 


WM»] 


Now,  let  E  be  of  the  form 


E  = 


(27) 


where  ek  is  defined  as 


= 


r  zrzlf  a2.  ....  «'"’■]  t=o . m- i 

C(Pk) 

P,*Pj 


(28) 


If  a  new  matrix  is  defined  as 
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C(p0)  0  •••  0 

Am  =  0  0  ,  (29) 

0  •••  0  C(pM_l)_ 

then  it  can  be  seen  that 

1  Po1  -  A(M)" 

i  i  i  i  ^ 

_1  p-J_,  p^!)_ 

i  a1  -  *rir 

=aaX  ;  :  :  :  .  oo) 

1  ox  ... 

_  HU- 1  HM- 1 

=  F 

r  MxM 

Substituting  the  results  of  (30)  into  (26),  we  reach  the  conclusion  that 

=  (31) 

The  transmitter-receiver  structure  in  (25)-(31)  is  more  general  than  standard 
OFDM.  It  can  be  easily  seen  that  the  choice 

'k2— 

pk=e~M  k  =  0,...,M  -1  (32) 

will  yield  the  standard  OFDM.  In  this  case,  the  operation  can  be  done  using  FFTs  if  the 
Rx  matrix  is  chosen  carefully;  namely, 

A,  =  VM\M  .  (33) 

2.  Channel  Estimation  from  a  Block  of  Received  Data 

If  the  received  vectors  are  organized  column-wise  into  a  matrix  and  the  number  of 
received  vectors  is  much  larger  than  the  length  of  the  OFDM  symbol  such  that 

Y  =  [y[0],  y[l],  ...,  y[J- 1]]  J»P,  (34) 

then  it  follows  that  Y  is  a  PxJ  matrix  of  the  form 

Y  =  ARt  [y[0],  5[1],  ...,  (35) 

From  (35)  it  can  be  seen  that  Y  is  not  a  full  rank  matrix  since  the  columns  of  A 
in  (22)  span  an  M  -dimensional  subspace  of  C1’ .  Therefore  ,a  matrix  V  exists  such  that 
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VY  =  0  (36) 

Since  the  matrix  of  transmitted  data  is  most  likely  full  rank,  the  matrix  V  also 
satisfies  VA-  0  or,  more  explicitly, 


V 


_  c[0] 

0 

0 

c[l] 

c[0]  ••• 

41]  - 

c[L~  1] 

0 

0 

c[L- 1]  ••• 

40] 

0 

41] 

0 

0 

c[L-l] 

=  0. 


(37) 


The  matrix  V ,  which  can  be  computed  from  the  null  space  of  Y ,  contains  the 
information  of  the  channel,  as  can  be  seen  next. 

To  extract  the  channel  information  from  the  V  matrix,  several  more  steps  are 
required.  Let  G  be  a  PxP  matrix,  featuring  a  structure  similar  to  E  in  (27), 


G  = 


(38) 


where 


=1 


42, 


pk 


(p- 1) 


? 


=  [!»  Pk1’  Pk 

k  =  0,...,P-l  .  (39) 

P  *  P} 

Although  this  is  more  general,  OFDM  would  correspond  to  the  choice  pk  =  e1"'1  for 
some  wk .  Now,  let  the  G  matrix  be  used  for  post-processing  of  the  received  signal.  For 
the  illustration  of  this  post-processing,  a  new  y'[n]  is  defined: 
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/[n\  =  GARls[n].  (40) 

Notice  that  (38)  and  (40)  yield  a  similar  receiver  to  the  one  defined  by  (25)  and 
(27),  only  with  an  extended  set  of  frequencies  at  the  demodulator. 


Now,  substituting  (38)  for  G  into  (40),  we  obtain 

T 

io 

T 


y'[n\  = 


~M-\ 

T 


(41) 


If  VLxM  is  defined  as  the  first  M  columns  of  G  ,  AM  is  defined  in  (29)  and  AL  is 


defined 


K  = 


then  (41)  can  be  written  as 


y\n\  = 


<(Pm) 

0 

0 

0 

0  c 

"Am 

0  " 

~v 

r  MxM 

0 

A,. 

v 

L  LxM  J 

Since  Rx  was  chosen  as  VMlxM  ,  (43)  can  now  be  rewritten  as 


(42) 


(43) 


y'M  = 


"A  M 

0  " 

0 

AI_ 

B 

y[//] . 


(44) 


where 


B  =  VLxMRx  (45) 

Recall  that  V  is  an  M  xP  complex  matrix  such  that  the  product  of  V  and  the  Y 
matrix  defined  in  (34)  is  equal  to  zero.  It  follows  that 

Vy[n\  =  0  V«  .  (46) 

This  implies 
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V 


"Am 

0  ' 

~I~ 

0 

Ai_ 

B 

=  0 


(47) 


or,  by  properly  partitioning  V , 


K  V2 ] 


Am 

A  lB_ 


=  0 


To  satisfy  (48),  the  V  matrix  can  be  chosen 

v  =  [K  v2]=[a lba~J  -/], 


(48) 


(49) 


where  I  is  an  L  x  L  identity  matrix.  This  V  matrix  must  also  be  equal  to  zero  when 
multiplied  with  the  Y  matrix,  meaning  that 


a^aa; 


^0 

Y, 


=  0. 


(50) 


where  Y0  is  a  matrix  of  the  first  M  rows  of  Y ,  and  Yx  is  a  matrix  of  the  last  L  rows  of 
Y  .  From  (50)  it  is  easy  to  see  that 


A  lBK~X  =  Yx. 


(51) 


Using  linear  algebra  on  (5 1),  we  get 

A  LBAlY0=Yl3 


alba-M  =  yX, 


T*  \  ^ 


Let 


KbK  =W'(W') 

hLB  =  Yfi{YX)-'  A 


(52) 


where 


Z„ 


nLxM 


Equation  (52)  then  becomes 


AlB  =  Z0Am- 


(53) 

(54) 

(55) 


If  the  structure  of  (55)  is  examined  more  closely  and  the  k  -th  row  of  both  sides 
are  taken, 
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¥  1  [c(a) 

C{Pm+]c  )  At,  0  “■  At.M- 1  =  [Z/t,0’  •••’  ZA,M-l]  ’•  ’(56) 

'■•I  J  L  0  c(pM- 1). 

then  \ 

•••’  At,M-l  ]  =  [^.O^A) )’  Zi,l^(A)’  •••’  •  (57) 

With  one  last  manipulation,  the  frequency  response  for  each  of  the  frequencies 
associated  with  the  M  point  IFFT  is 

C(A)  = - —  Li 

C(pMtl)zu 

l  =  0, 1,...,M  -1  .  (58) 

k  =  0,1,...,L-1 

The  C(pM+k )  term  is  a  constant  associated  with  all  frequencies  of  the  channel 

response.  Additionally,  for  robustness  in  the  presence  of  noise,  the  channel  response  can 
be  averaged  over  the  k  rows  of  the  matrix. 

In  conclusion,  the  channel  frequency  response,  apart  from  a  scaling  factor,  can  be 
computed  from  a  block  of  received  OFDM  symbols,  as  in  (34),  as  follows: 

Step  1 :  Let  VLxM ,  VMxM  be  the  first  LxM  and  the  last  LxL  blocks  of  G  in 
Equation  (38). 

Step  2:  Computed  =  VLxMRx . 

Step  3:  Let  Y0  and  Yx  be  the  first  M  rows  and  last  L  rows  of  Y  ,  respectively. 

Step  4:  Compute  Z0  =  YXY*  ( Y(l Y*  j  . 

Step  5:  Call  zk  and  bk  any  row  k  of  Z0  and  B  . 

Step  6:  C(p0),  C(px),  ...,  C{pM_x)  is  proportional  to  ratio  of  rows  zk ^  and  bh .. 

The  channel  frequency  response  was  calculated  using  a  MATLAB  program  that 
implemented  the  above  method.  The  next  chapter  discusses  the  adaptive  method  used  for 
real-time  processing  of  an  OFDM  signal. 
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IV.  RECURSIVE  IMPLEMENTATION 


A.  OFDM  WITH  CYCLIC  PREFIX  AND  OFDM  WITH  ZERO  PREFIX 

The  blind  channel  equalization  algorithm  presented  in  Chapter  III  is  not  recursive, 
and  a  large  number  of  received  OFDM  symbols  are  required  to  compute  the  channel 
estimate.  In  this  chapter,  a  recursive  method  of  computing  the  channel  estimate  for  every 
block  received  using  an  LMS  estimator  is  explained  in  detail.  For  the  purposes  of  this 
portion  of  the  thesis,  a  received  signal  similar  in  structure  to  that  used  by  IEEE  802. 1  la  is 
used.  In  this  standard,  the  length  of  the  IFFT  M  is  64  samples,  and  the  length  L  of  the 
cyclic  prefix  is  16,  for  an  overall  OFDM  symbol  length  P  of  80  samples  [6]. 

The  method  of  blind  equalization  described  in  Chapter  III  required  a  ZP.  In  this 
chapter,  we  will  see  that  the  CP  OFDM  symbol  used  in  all  standards  can  be  reframed  as  a 
ZP  OFDM  symbol  by  simple  preprocessing  at  the  receiver. 

To  achieve  this  goal,  a  standard  CP  OFDM  signal,  yCP[n\ ,  can  be  made  into  a  ZP 
OFDM  signal  by  subtracting  a  delayed  version  of  the  signal  to  get 

yzp\.1A  =  yCp\.nl\-ycp\.n-M1\  ■  (59) 

As  shown  in  Figure  4,  it  is  easy  to  see  that  the  CP  terms  cancel  since  the  data  in 
the  CP  is  the  same  as  the  data  at  the  end  of  the  OFDM  symbol.  To  recover  the  old 
symbols,  the  new  OFDM  symbol  and  the  phase  shifted  old  symbol  are  simply  added. 


CP  1 

OFDM  SYMBOL  1 

CP2 

OFDM  SYMBOL  2 

CP  1 


OFDM  SYMBOL  1 


CP  2 


OFDM  SYMBOL  1  WITH 
PHASE  SHIFT 


ZP 


NEW  OFDM  SYMBOL 


Figure  4.  A  delay  of  M  points  results  in  a  phase  shift  of  the  first  OFDM  symbol  with 

a  ZP  and  new  symbols  subsequent  to  that. 
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Now  the  channel  can  be  estimated  using  the  algorithm  discussed  in  the  previous 
chapter  and  applied  to  the  received  OFDM  symbol.  The  rest  of  this  section  describes  the 
LMS  filter  method  used  as  if  a  zero  prefix  OFDM  symbol  had  been  transmitted  vice  the 
CP  OFDM  symbol  that  was  actually  simulated. 

Initially,  an  LMS  estimator  was  derived  directly  from  the  off-line  method 
described  in  Chapter  III.  This  method  is  not  robust,  and  the  theory  behind  it  is  not 
presented  in  this  thesis;  however,  some  of  the  simulation  results  have  been  included  for 
comparison  purposes.  In  what  follows  we  present  a  different  approach  to  blind 
equalization  for  ZP  OFDM,  discussed  in  [9],  which  is  more  suitable  to  recursive 
implementation. 


B.  PROPOSED  APPROACH 


In  this  section,  a  different  approach  to  blind  equalization  of  ZP  OFDM  is 
investigated  based  on  standard  properties  of  the  DFT.  The  results  of  this  section  are 
based  on  [9]. 

The  purpose  of  this  approach  is  an  improved  estimation  method  that  can  be 
computed  recursively  with  an  LMS.  The  resulting  method  is  shown  to  be  effective  in 
cases  where  the  channel  is  initially  poorly  estimated  or  slowly  changing  with  time.  It 
also  had  an  added  benefit  that  is  discussed  in  the  next  section  of  this  chapter. 


In  order  to  introduce  this  approach,  define  a  vector  x  of  length  2 M  such  that 

x[n\  =  0,  for n  =  M,...,2M (60) 
and  a  second  vector  w  of  the  same  length  2 M  of  the  form 

f  0  if  0  <  n  <  M  - 1 

u{»]=  •  (61) 

[l  if  M  <n<  2M  - 1 

Clearly,  no  matter  what  x[n\,  n  =  0 -1  is,  we  have 

x[ri\w[n\  =  0  for  all  n  =  0,...,2M -1 .  (62) 

Since  multiplication  of  two  finite  length  sequences  corresponds  to  the  circular 
convolution  of  the  respective  DFTs,  we  can  write  (62)  in  terms  of  W[k]  =  DFT {w[/z]] 

and  X[k]  =  DFT {x[n§ . 
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The  DFT  of  the  w  vector  is  defined  as 


2M-\ 

W[k ]  =Yae  1M  ’  (63) 

n=M 

which  is  a  geometric  series.  If  l  =  n-M^>n  =  l  +  M,  then  (63)  can  be  rewritten  as 

for  k  =  -1 


W[k]  = 


V-I  -j—m+M)  ,  l-(-l) 

Y^e  j2n  =(_!)*  1 


1= 0 

M,  for  k  =  0 


-j—k 

l-e  N 


(64) 


It  is  easy  to  see  that  W[k\  =  0  when  k  is  even,  and  in  particular 

-2 


W[k]  = 


-j—k 

l-e  N 


0 

M 


if  k  odd 

if  k  *  0,  even . 
if  A:  =  0 


(65) 


Now,  circularly  convolve  the  DFTs  of  the  two  vectors  x  and  w  to  get 

M- 1 

=  0,  (66) 

1=0 

where  the  shifted  represented  by  k-l  is  a  circular  shift.  If  k  is  chosen  to  be  any  odd 
frequency,  it  can  be  shown  that  W[k  - 1 ]  is  not  equal  to  zero  only  when  k  =  l  or  /  is  even. 
Using  this  conclusion,  we  see  that  when  k  is  odd  (66)  becomes 

M- 1 

Jl[/F[U/]  +  I[iM0]  =  0.  (67) 

1=0 
l  even 

If  the  value  for  W[0]  is  substituted  into  (67)  and  the  equation  is  rearranged,  then 

M- 1 

MX[k]  =  ~Y,  W[k  -  l]X[l\ .  (68) 

/=0 
l  even 

The  relationship  between  the  DFT  of  the  transmitted  signal  and  the  DFT  of  the 
received  signal  is  just  the  inverse  of  the  channel  frequency  response;  that  is, 

X[k\  =  ^-.  (69) 

C[k\ 

For  any  odd  value  of  k ,  (69)  can  then  be  substituted  into  (68)  to  get 
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MY[k] 


C[k] 


M  — 1  1 

Y'w-mn— 


1=0 

y  /  eve« 


(70) 


If  (70)  is  further  deconstructed,  with  the  /  values  associated  with  the  data  on  the 
left-hand  side,  the  /  values  associated  with  the  pilots  on  the  right-hand  side,  and  an 
arbitrary  odd  value  of  k ,  then 


M-l  1  1 

Y  W[k-l]Y[l] - +  MY[k ] - 

^  C[l]  C[k] 


1=0 
l  even 


(  \ 
M- 1 

1=0 

y  /  even 


(71) 


With  a  few  terms  further  defined,  (71)  lends  itself  to  an  LMS-type 
implementation;  that  is, 

£n=[W[k]Y[  0],  W[(k-2)]Y[2\,  ...,  MY[k]] 

11  11 


6  = 


C[0]  C[2] 

c 


C[M-  2]  C[A:] 


(72) 


y^-/]xp,oj/] 


/=o 

y  /  even 


Now,  from  [9],  the  LMS  estimator  can  be  written  as 


0  =  0  i .  +  Lid  ( 0Tlr, (h  - v, 

— new  — old  *Z_k  n  \ — °'“—k  n  k,n  J 


(73) 


C.  SIMULATION  RESULTS 

Each  of  the  methods  of  implementation  for  the  blind  estimation  of  the  channel 
were  run  using  100  trials  at  each  data  point,  regardless  of  the  type  of  test  being  run. 
Additionally,  the  number  of  OFDM  symbols  transmitted  was  set  at  10,000.  Using  IEEE 
Standard  802.1  la  format,  we  have  52  BPSK  or  QPSK  symbols  per  OFDM  symbol,  for  a 
total  of  52  million  symbols  per  data  point. 

The  number  of  pilot  signals  was  reduced  from  four  pilot  symbols  per  OFDM 
symbol  to  one  pilot  symbol  per  five  OFDM  symbols.  The  purpose  of  transmitting  the 
pilot  symbol  was  to  determine  the  channel  constant  for  all  three  methods  of 
implementation.  However,  the  number  of  pilot  signals  sent  could  be  further  reduced 
depending  on  channel  conditions. 


24 


The  values  of  // ,  the  LMS  algorithm  step  size,  were  determined  experimentally  to 
ensure  convergence  of  the  adaptive  algorithm.  For  the  initial  LMS  implementation  the 
step  size  //  was  set  to  a  small  value  ( 1 0  7 ),  and  for  the  second  DFT-based  LMS 
implementation,  //  was  set  to  a  larger  value  (3x1 0  5 ).  These  values  can  be  adjusted  to 
ensure  the  algorithm  converges  or  to  adjust  the  speed  at  which  the  algorithm  converges. 

1.  Performance  in  the  Presence  of  Multipath  and  AWGN 

To  test  the  performance  of  the  proposed  DFT-based  adaptive  channel  estimation 
method  and  the  adaptive  method  based  on  the  theory  discussed  in  Chapter  III,  the 
respective  performances  were  compared  using  the  initial  off-line,  or  block,  blind 
estimation  method.  All  three  methods  were  tested  with  AWGN.  The  tests  were  run  with 
SNRs  from  zero  to  50  dB.  Algorithm  perfonnance  was  determined  by  calculating  the 
symbol  error  rate  (SER),  in  errors  per  symbol.  The  results  for  the  DFT  based  adaptive 
algorithm  with  an  LMS  step  size  of  //  =  1 0  and  the  block  method  is  shown  in  Figure  5. 
The  channels  are  modeled  as  random  complex  channels  with  an  impulse  response  of  the 
same  length  as  the  guard  interval  (L  =16  in  this  case). 
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SNR(dB) 


Figure  5.  QPSK  SER  versus  SNR  with  a  step  size  of  //  =  10  7 . 


For  this  thesis  the  SNR  is  defined  as 


SNR  = 


E  [s\t)] 
E  Wit)} 


where  s(t )  is  the  signal  and  n(t)  it  the  noise.  It  is  related  to  Eh/Na 
where  B  is  the  bandwidth  of  the  channel,  and  Rh  is  the  data  rate, 
this  factor  is  0.692  [6]. 


(74) 

by  a  factor  of  B/Rb 
For  802.11a  QPSK, 


It  can  be  seen  that,  in  this  example,  the  proposed  LMS  method  has  better 
performance  than  the  block  processing  approach  presented  in  Chapter  III  for  SNR  less 
than  18  dB.  At  5  dB,  the  proposed  method  out  performs  the  block  method  by 
2.423x10  1  errors/symbol,  the  maximum  delta  between  the  two  methods.  This  is  likely 
due  to  the  fact  that  the  block  processing  approach  estimate  is  averaged  over  the  k  rows 
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of  the  matrix,  as  discussed  in  the  previous  chapter.  Specifically,  there  are  k  estimates  of 
the  channel,  one  derived  from  every  row  of  the  matrix,  and  these  estimates  are  averaged 
together  before  being  applied  to  the  received  signal. 

The  results  for  the  adaptive  method  discussed  in  Chapter  III  are  not  presented 
since  the  algorithm  never  fully  converged  to  the  correct  channel.  This  is  due  to  the  fact 
that  the  matrices  that  the  LMS  method  is  based  on  are  ill-conditioned.  The  block  method 
uses  singular  value  decomposition  to  overcome  this  limitation. 

The  next  simulation  was  done  under  the  same  conditions  except  the  value  of  the 
step  size  was  changed  to  /u  =  3x10  ''  to  determine  the  effect  on  the  proposed  algorithm. 
The  results  of  this  change  can  be  seen  in  Figure  6.  The  proposed  adaptive  method  works 
poorly  compared  to  the  block  estimation  method. 


Figure  6.  QPSK  SER  versus  SNR  with  a  step  size  of  ju  =  3x10  5 . 
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2.  Non-Stationary  Channel  Performance 

The  next  test  of  the  proposed  implementation  was  to  examine  the  effect  that  a 
non-stationary  channel  has  on  its  performance.  The  non-stationary  channel  was 
simulated  using  a  basic  random  walk.  The  channel  was  initially  a  set  of  four  complex 
impulses  distributed  within  the  length  of  the  guard  interval.  Every  OFDM  symbol,  this 
channel  was  updated  to  a  new  value  by  a  random  walk.  The  random  walk  step  size  has 
been  calculated  as  a  percent  of  the  channel  power.  The  channel  coefficients  are  obtained 
iteratively  from 

Ct+l  =ct+  /?E[C,2]AC  (75) 

where  AC  is  a  random  vector  consisting  of  [l,  0,  — l] ,  /?  is  a  scaling  factor,  and  E[C,2]  is 
the  channel  power. 

The  block  method  of  channel  estimation  was  expected  to  perform  poorly,  based 
on  the  fact  that  it  assumes  a  stationary  channel  over  the  entire  data  set  of  10000  received 
blocks.  With  the  small  LMS  step  size  ju  =  1 0  7 ,  the  adaptive  method  performed 
comparable  to  the  block  method.  The  small  step  size  affected  its  ability  to  converge  on 
the  actual  channel  estimate  as  each  OFDM  symbol  was  received. 

To  ensure  that  the  poor  perfonnance  was  due  to  the  LMS  step  size,  the  simulation 
was  also  run  using  //  =  3  x  1 0  5 .  The  results  of  this  simulation  are  shown  in  Figure  7, 
where  SNR  is  defined  in  Section  1  of  this  chapter. 
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Figure  7.  QPSK  SER  versus  random  walk  step  size  with  an  LMS  step  size  of 

ju  =  3x  1CT5 . 

As  illustrated  in  Figure  7,  the  LMS  method  based  on  the  blind  equalization 
algorithm  in  Chapter  III  is  very  sensitive  to  the  LMS  step  size.  The  proposed  DFT  based 
adaptive  method  performs  better  than  the  block  method  and  the  Chapter  III  LMS  method 
fails  to  converge  to  a  solution.  The  step  size  for  the  proposed  DFT  based  adaptive 
method  can  be  adjusted  to  give  faster  convergence  on  a  non-stationary  channel;  however, 
the  gains  are  marginal,  since  the  block  method  reaches  10%  SER  for  a  random  walk  step 
size  of  0.21  percent  of  channel  power,  while  the  alternative  LMS  method  reaches  10% 
SER  at  0.23  percent.  This  can  also  be  seen  in  Figure  8,  where  the  alternative  LMS 
method  and  block  method  have  been  plotted  using  a  logarithmic  scale.  The  error  rate  for 
the  block  method  is  better  than  that  of  the  alternative  LMS  method  until  a  random  walk 
step  size  of  0.08  percent  of  the  total  channel  power  is  reached.  Although  there  are  several 
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data  points  at  which  the  block  method  outperforms  the  alternative  LMS  method  above 
0.08  percent,  in  general,  the  alternative  LMS  performs  better. 


Figure  8.  A  comparison  of  the  block  and  proposed  adaptive  method  of  calculating  the 
channel  using  a  step  size  of  //  =  3  x  10  5 .  These  results  were  generated 
using  500  trials  at  each  data  point. 


These  results  also  fail  to  take  into  account  the  fact  that  the  larger  step  size  affects 
the  performance  in  AWGN,  as  seen  in  Section  1. 

3.  Effects  of  Errors  in  the  Initial  Channel  Estimate 

The  simulations  performed  in  the  preceding  sections  assumed  that  there  was  a 
perfect  initial  estimate  of  the  channel  conditions  based  on  the  reception  of  a  preamble. 
The  performances  of  the  algorithms  based  on  a  poor  initial  channel  estimate  are 
investigated  in  this  section.  Only  the  two  LMS-based  methods  are  tested  in  this  section 
since  the  block  method  does  not  rely  on  an  initial  channel  estimate.  The  error  in  the 
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initial  estimate  was  simulated  by  taking  the  original  channel  estimate  and  multiplying 
each  frequency  of  the  channel  by  a  different  random  real  number.  The  initial  estimate  is 
given  by: 

CinitJ  =  CacjVj  (76) 

where  Cinit  is  the  channel  estimate  given  to  the  simulation,  Cact  is  the  actual  channel 


impulse  response,  and  v.  ~  /V(l,r)  with  r  defined  as  the  seed  channel  variance.  The 

performance  metric  for  this  algorithm  was  based  on  the  number  of  OFDM  symbols  it 
took  for  the  algorithm  to  converge  to  less  than  a  10%  SER.  This  was  calculated  by 
running  the  algorithm  until  an  OFDM  symbol  had  five  or  less  errors  in  it.  The  number  of 
symbols  was  still  set  to  10,000,  and  if  the  algorithm  failed  to  converge  by  then,  it  was 
assumed  that  it  would  not  converge.  The  results  of  this  simulation  are  seen  in  Figure  9. 


Figure  9.  Number  of  OFDM  symbols  required  to  converge  versus  seed  channel 

variance  for  a  step  size  of  /z  =  1 0  7 . 
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Both  algorithms  perform  poorly  at  the  lower  step  size  of// =  10  7 ,  taking  over 
8,000  blocks  to  converge  at  a  seed  channel  variance  of  1.0  and  failing  to  converge 
beyond  a  seed  channel  variance  of  1.4.  Assuming  a  OFDM  symbol  interval  of  4.0  jus, 
based  on  IEEE  Standard  802.11a,  we  find  that  this  translates  to  a  convergence  time  of  32 
ms. 

The  step  size  was  then  increased  to  //  =  3xl0~5  to  see  what  effect,  if  any,  it  had 
on  the  convergence  of  the  two  algorithms.  It  can  be  seen  in  Figure  9  that  the  alternative 
LMS  method  converges  much  faster  than  the  original  LMS  method. 
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Figure  10.  Number  of  OFDM  symbols  required  to  converge  versus  seed  channel 

variance  for  a  step  size  of  //  =  3 xlO  5 . 
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With  the  larger  step  size,  the  alternative  LMS  method  converges  between  20  and 
24  ms  for  seed  channel  variances  of  1.1  to  2.0.  The  original  LMS  method  shows  slightly 
degraded  perfonnance  at  the  larger  step  size.  Due  to  the  fact  that  the  algorithms  do  not 
always  converge,  the  plot  is  not  smooth  through  seed  channel  variances  of  1.3  to  2.0. 

This  chapter  investigated  a  DFT-based  adaptive  method  for  channel  estimation. 
The  next  chapter  discusses  an  algorithm  for  the  estimation  of  data  in  null  subchannels 
that  is  also  based  on  the  DFT. 
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V.  DATA  RECOVERY  FROM  A  FADED  SUBCARRIER 


A.  PROPOSED  NULL  ESTIMATION  METHOD 

One  of  the  issues  in  OFDM  is  the  loss  of  data  in  subcarriers  corresponding  to 
nulls  in  the  frequency  response  of  the  channel.  This  is  particularly  significant  in 
environments  with  multipath,  which  adds  destructively  at  certain  frequencies.  The 
method  presented  in  this  chapter  is  based  on  [10]. 

The  argument  of  this  section  is  that,  using  OFDM  with  ZP  and  an  FFT  size  twice 
the  OFDM  symbol  length,  the  even  and  odd  subcarriers  of  the  OFDM  symbol  offer 
sufficient  redundancy  so  that  the  lost  data  can  be  recovered. 

In  order  to  introduce  this  concept,  let  the  transmitted  OFDM  symbol  be 

x  =  [v[0],v[l],...,x[M-l],0,...,0]  (77) 

where  x  is  a  vector  of  length  P ,  which  includes  the  ZP  of  L  zeros  at  the  end.  Now 
define  the  M  point  DFT  and  the  2 M  point  (zero  padded)  DFT  as 
XM[k]  =  DFT{[x[0],...,x[M -l]]},  k  =  -1 

X2M [m\  =  DFT{[x[0],...,x[M  - 1]] , 0, 0, . . . , 0} ,  m  =  0,...,2M -1  ’ 
where  M  zeros  are  added  to  x  before  the  2 M  point  DFT  is  taken. 


Then,  as  discussed  in  (66),  for  any  odd  k  we  can  write 

2M-\ 

£  X2M  [l]W[k  -  /]  +  X2M  [k]W[ 0]  =  0 .  (79) 

1=0 
l  even 

Since  t  is  even  and  k  is  odd  we  can  let  k  -  2m  + 1  and  rewrite  this  equation  as 

M  JP”* + ']  =  -Y,  W[2 (m  - 1)  + 1  ]XM  [f] .  (80) 

C[2m  + 1]  to 

Now,  set  an  arbitrary  threshold  and  use  it  to  partition  the  transmitted  data  into  two 
parts  such  that: 

X[l]  =  x+[l]  +  x_[l]  (81) 


with 
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\X[l]  if  \C[l]\>  Thresh 
0  if  \C[l]\<  Thresh 

\X[1 ]  if  \C[I]\<  Thresh 
0  if  \C[I]\>  Thresh 

where  Thresh  is  an  arbitrary  threshold. 


X+[f  = 
Xfl]  = 


(82) 


In  other  words,  X_[f,]  indicates  the  subcarriers  received  at  low  power  (at  or  close 
to  a  null  in  the  channel)  and  X+[f]  all  others. 


Now,  if  the  partitioned  data  is  inserted  into  (80),  it  becomes 
Y\  + 11  1  M  1  M  1 

- 77  =  —  Z  w - 0  +  1R [/]  +  —  X  w^m  -0  +  IK [/]  • 

C[2m  + 1]  Mto  Mto 


(83) 


-J — n 
M 


(84) 


Define 

w0[n]  =  IDFT{W  [2/c  +  l]}  =  -e 
x ±[n\  =  IDFT{X±[l]} 

Using  the  above  definitions  in  (83),  we  get 

^2k++X\  =  DFT{wMxM}+DFT{wfn]x_[n]}  . 

If  the  term  containing  x_[»]  on  the  right  hand  side  of  the  equation  is  isolated,  then 


(85) 


DFT{wfn\x_[n]}  =  - DFT {w0[n]x+[n]}  . 


(86) 


C[2k  + 1] 

If  the  inverse  discrete  Fourier  transform  (IDFT)  of  both  sides  of  (86)  is  taken 


then 


wfn\x_\n\  =  IDFT 


Y\2k  +  1] 


-w0[n]x+[n\. 


(87) 


[C[2k  +  1] 

Now,  dividing  both  sides  by  the  w0[n]  tenn  to  solve  for  the  terms  of  the 
transmitted  signal  below  the  threshold,  we  get 


and,  finally, 


x  [//]  =  — —  IDFT  {  7[2A  +1]1  -  x  [/?] 
-L  w0[n]  \  C[2k  + 1]  J 
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(88) 


(89) 


X_[e]  =  DFT{x_[n]},l  =  0,...,M-l 
is  the  data  recovered  at,  or  close  to,  the  nulls  of  the  channel. 

B.  SIMULATION  RESULTS 

To  test  the  null  estimation  method  described  in  Section  A  of  this  chapter,  a  test 
similar  to  those  in  Chapter  IV  were  run,  with  several  significant  differences.  Instead  of 
using  the  OFDM  symbol  fonnat  described  in  IEEE  Standard  802.11a,  no  nulls  or  pilots 
were  included  in  the  OFDM  symbol.  This  was  done  in  order  to  test  the  proposed  concept 
of  null  estimation  without  having  to  worry  about  the  effects  of  the  nulls  in  the  OFDM 
symbol.  Additionally,  knowledge  of  the  channel  was  assumed. 

The  null  estimation  algorithm  continued  to  use  the  delay  of  one  block  to  achieve 
the  effects  of  a  ZP  while  transmitting  the  symbol  with  a  CP.  To  test  the  algorithm, 
transmission  was  simulated  over  a  noisy  channel,  where  the  noise  was  AWGN. 

For  testing,  we  used  a  family  of  random  channels  all  with  zero  frequency  response 
at  a  specific  frequency  <x>0 .  An  example  channel  response  is  shown  in  Figure  11. 
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Figure  11.  One  of  a  family  of  channel  frequency  responses,  all  having  nulls  at 

co  =  0.589  radians. 
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Figure  12.  A  comparison  of  the  QPSK  symbol  error  rate  between  a  normal  OFDM 
receiver  algorithm  and  a  null  estimating  OFDM  receiver  algorithm. 

The  performance  of  the  standard  algorithm  (solid  line)  is  compared  with  that  of 
the  proposed  algorithm  (dashed  line)  in  terms  of  symbol  error  versus  SNR  in  Figure  12. 
The  two  algorithms  have  similar  performance  up  to  25  dB  SNR,  although  the  null 
estimating  OFDM  receiver  algorithm  performs  slightly  better  overall.  Above  25  dB 
SNR,  there  is  a  sharp  drop  in  the  error  rate  from  the  null  estimation  algorithm,  while  the 
normal  OFDM  receiver  algorithm  remains  constant  at  7.8x10  errors  per  symbol  due  to 
the  information  lost  in  the  deep  fading  sub-carriers. 

This  chapter  presented  an  algorithm  for  the  estimation  of  the  data  lost  in  null 
subchannels  and  simulation  results  from  the  implementation  of  this  algorithm.  Chapter 
VI  provides  a  summary  of  the  work  presented  in  this  thesis,  significant  results,  and 
recommendation  for  future  work. 


39 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


40 


VI.  CONCLUSION 


The  goal  of  this  thesis  was  to  investigate  an  adaptive  implementation  of  the  blind 
channel  equalization  method  proposed  in  [1],  This  adaptive  method  was  implemented 
using  the  current  OFDM  symbol  architecture  from  IEEE  standard  802.11a.  Once  the 
adaptive  method  based  on  the  work  in  [1]  was  simulated,  an  alternative  method  was 
developed  with  more  robust  perfonnance  characteristics.  The  algorithm  for  the  second 
method  led  to  the  development  of  a  proposed  method  of  null  channel  estimation,  which 
was  also  simulated  in  this  thesis. 

A.  SUMMARY 

This  thesis  presented  a  subspace-based  method  of  blind  channel  estimation  of 
OFDM  modulated  signals  that  used  the  zero  prefix  as  the  guard  band  interval.  Based  on 
this,  a  method  of  adaptive  blind  channel  estimation  has  been  proposed.  Furthermore,  it  is 
shown  that  the  redundancy  in  ZP  OFDM  can  be  exploited  to  recover  data  from  deeply 
faded  subcarriers.  The  blind  channel  estimation  methods  were  implemented  in  MATLAB 
and  then  compared  for  perfonnance  for  several  different  channel  conditions.  Both  a 
noisy  channel  with  AWGN  and  a  non- stationary  channel  were  used  to  compare  the 
results  of  the  three  methods  of  channel  estimation.  Additionally,  two  adaptive  methods 
were  tested  in  the  presence  of  enors  in  the  initial  channel  estimate  to  see  if  they  could 
recover  the  original  channel  estimate. 

Finally,  the  algorithm  for  data  recovery  from  faded  subcarriers  was  coded  in 
MATLAB  and  simulated.  Since  in  this  case  we  assume  full  knowledge  of  the  channel,  it 
was  only  tested  in  a  channel  with  AWGN.  The  results  of  this  simulation  were  compared 
to  those  of  an  OFDM  system  that  did  not  perform  null  estimation. 

B.  SIGNIFICANT  RESULTS 

The  method  of  blind  channel  estimation  proposed  in  [1]  can  be  implemented  in  an 
adaptive  fonn.  However,  this  method  is  very  sensitive  to  the  step  size  used  in  the  LMS 
algorithm.  The  algorithm  did  not  converge  in  a  channel  with  AWGN  with  a  relatively 
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large  step  size  (//  =  3xlO  5  in  the  example).  Similarly,  the  algorithm  performed  poorly 
with  this  step  size  for  both  the  non-stationary  channel  simulation  and  the  incorrect  initial 
channel  estimate  simulation. 

The  proposed  adaptive  method  is  more  robust  than  the  original.  For  the  AWGN 
noisy  channel  with  an  LMS  step  size  of  //  =  1 0  7 ,  the  system  converged  to  a  lower  error 
rate  in  simulated  low  SNR  environments  than  the  off-line,  block  method.  This  can  be 
seen  in  Table  1. 


Table  1.  Comparison  of  calculated  SER  at  various  SNR  for  the  blind  estimation 
methods  using  a  step  size  of  /z  =  1 0  7 . 


SNR  (dB) 

Block  Method 
(Errors/Symbol) 

Proposed  DFT 
Based  Adaptive 
Method 

(Errors/Symbol) 

0 

5.799x10"' 

4.543x10"' 

5 

4.513x10"' 

2.090x10"' 

10 

1.741x10  1 

6.203  xlO"2 

15 

3.895xl0"2 

2.155xl0’2 

20 

2.640x1 0'3 

5.982  xlO"3 

25 

1.390xl0"5 

4.837  xlO’4 

With  the  larger  step  size,  the  block  method  outperforms  the  adaptive  LMS 
algorithm.  However,  there  are  several  advantages  to  the  larger  LMS  step  size,  which  can 
be  seen  in  the  results  from  the  other  two  simulations.  The  alternative  adaptive  method 
with  the  larger  step  size  outperfonns  the  block  algorithm  during  the  non-stationary 
channel  simulation  above  a  random  walk  step  size  of  0.08  percent  of  channel  power  by  an 
average  of  0.0175  symbol  errors/symbol.  For  the  bad  initial  channel  estimation 
simulation,  the  larger  LMS  step  size  outperfonns  the  smaller  step  size  significantly, 
which  can  be  seen  in  Table  2. 
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Table  2.  Comparison  of  convergence  speed  for  various  seed  channel  variances  for 
the  original  and  alternative  adaptive  methods  with  an  LMS  step  size  of  //  =  3x  1CT5 


Seed  Channel 
Variance 

Adaptive  Method 
(Number  of  OFDM 
Symbols) 

Alternative  Adaptive  Method 
(Number  of  OFDM 
Symbols) 

Time 

Difference 

(ms) 

0 

1 

1 

0 

0.5 

1001 

4.25 

4 

1.0 

9800 

4009 

23.2 

1.5 

Fails  to  Converge 

5049 

00 

2.0 

Fails  to  Converge 

6031 

00 

Perhaps  the  most  significant  result  is  the  fact  that  the  proposed  adaptive  method 
can  be  tuned  for  performance  based  on  the  channel  conditions  that  may  exist. 

The  recovery  of  data  transmitted  over  deep  fading  channels  was  shown  to  be 
possible  using  the  algorithm  described  in  Chapter  V.  For  relatively  high  SNR  (25  dB  in 
this  example),  the  proposed  algorithm  clearly  outperformed  the  standard  receiver  which 
cannot  recover  the  faded  data. 

C.  RECOMMENDATIONS  FOR  FUTURE  WORK 

There  are  numerous  areas  available  for  future  work  with  this  thesis.  First,  the 
proposed  algorithm  should  be  evaluated  using  actual  channel  characteristics  to  determine 
if  it  is  viable  for  use  in  an  actual  communication  method. 

Secondly,  the  adaptive  method  should  be  tuned  to  perform  based  on  channel 
characteristics.  In  a  nonnal  OFDM  modulation  scheme,  such  as  IEEE  Standard  802.11a, 
a  known  preamble  is  sent  before  data  transmission  begins.  The  channel  characteristics 
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determined  from  the  preamble  should  then  be  used  to  determine  the  desired  LMS  step 
size.  By  modifying  the  step  size,  the  performance  of  the  algorithm  can  be  enhanced  for 
noisy  or  non-stationary  channels. 

Third,  the  null  estimation  algorithm  should  be  tested  using  an  OFDM  symbol 
from  a  current  standard  and  integrated  with  channel  estimation  method.  If  knowledge  of 
the  channel  were  known  before  transmission,  no  data  would  be  sent  over  channels  with 
deep  fading.  Once  a  channel  has  been  determined  to  have  sub-carriers  with  deep  fading 
characteristics,  null  estimation  should  be  perfonned  to  determine  if  the  algorithm  will 
perform  as  well  with  only  an  estimation  of  the  channel  characteristics. 
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